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The effects of mode-mode and isospin-isospin correlations on nonequilibrium chiral dynamics are 
investigated by using the method of the time dependent variational approach with squeezed states as 
trial states. Our numerical simulations show that large domains of the disoriented chiral condensate 
(DCC) are formed due to the combined effect of the mode-mode and isospin-isospin correlations. 
Moreover, it is found that, when the mode-mode correlation is included, the DCC domain formation 
is accompanied by the amplification of the quantum fluctuation, which implies the squeezing of the 
state. However, neither the DCC domain formation nor the amplification of the quantum fluctuation 
is observed if only the isospin-isospin correlation is included. This suggests that the mode-mode 
coupling plays a key role in the DCC domain formation. 

PACS numbers: 25.75.-q, 11.30.Rd, ll.10.Lm 



I. INTRODUCTION 

The possibility of the formation of the disoriented chi- 
ral condensate (DCC) in relativistic heavy ion collisions 
has been investigated in a number of theoretical studies. 
Non-equilibrium field dynamics has been studied also in 
conjunction with the structure formation in early uni- 
verse. Within the classical approximation it has been 
recognized that large domains of DCC are produced dur- 
ing the nonequilibrium process in the course of the time 
evolution from the quench initial condition jj, . On the 
other hand, in calculations including quantum mechan- 
ical effects, mainly homogeneous (translationally invari- 
ant) systems have been studied to avoid technical diffi- 
culty in numerically calculating the two-point functions 
0, Q. In fact, there have been several attempts to in- 
clude spatial inhomogeneity H, IE an d also memory 
effect d, in the quantum mechanical time evolution 
of the fields. However, the correlation between modes 
with different momenta (mode-mode correlation) was not 
taken into account in these works and the effect of inter- 
actions has been included only through the mean fields. 
It has not been settled whether large domain structure is 
formed when quantum effects are taken into account. 

Recently, it was explicitly shown that DCC domains 
are formed in quantum calculation for the first time in the 
case of 1+1 dimensional geometry [Io|. In that calcula- 
tion, the mode-mode correlations are explicitly taken into 
account and no translational invariance was assumed. 
The importance of the mode-mode correlation is qual- 
itatively understood as follows. If one does not include 
the mode-mode correlation, each mode is decoupled from 
each other during the time evolution except for the in- 
direct interactions through the mean fields. The energy 
transfer from modes with short wavelengths to those with 
long wavelengths is not effective enough to form large cor- 
related domains. If one includes the mode- mode correla- 



tion, direct energy transfer from short wavelength modes 
to long wavelength modes and as a result the amplifica- 
tion of long wavelength modes in the pion fields become 
possible, which can lead to the formation of large DCC 
domains. 

In this paper, in addition to the mode-mode cor- 
relation, we examine the effect of the correlation be- 
tween fields having different isospin components (isospin- 
isospin correlation). This effect on the domain formation 
of the DCC has not been given particular attention to in 
earlier studies. However, the isospin-isospin correlation 
can arise during the nonequilibrium process as well as 
the mode-mode correlation. It is important to incorpo- 
rate both correlations on an equal footing. In addition, 
it is of interest to study whether the isospin-isospin cor- 
relation affects the domain formation of DCC and, if yes, 
it is then important to compare the ways the two correla- 
tions affect the domain formation. We also make further 
examination on the quantum mechanical features of the 
domain formation in this paper. 

The rest of the paper is organized as follows. In Section 
[n]we briefly review the formalism and the equations of 
motion, and discuss the initial condition. In Section IlIII 
we present numerical results. Section I1VI is devoted to 
summary and conclusions. 



II. EQUATIONS OF MOTION AND INITIAL 
CONDITION 

We take the 0(4) linear sigma model as a low energy ef- 
fective theory of QCD and apply the method of the time 
dependent variational approach (TDVA) with sq ueezed 
states as trial states as in the previous paper 10]. This 
method was originally developed by Jackiw and Kerman 
as an approximation in the functional Schrddinger ap- 
proach [Tl| and later it was shown to be equivalent to 
TDVA with squeezed states by Tsue and Fujiwara [T^ |. 
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This method makes it possible to solve the time-evolution fluctuations and correlations in a self-consistent manner, 
of the order parameters (mean fields), and the quantum 

I 



A. Equations of motion 



We denote the sigma (a(x)) and pion (n(x)) field operators as a four dimensional vector, <p a (x) — (a(x),T?(x)), 
where a runs from to 3. Then the Hamiltonian H of the 0(4) linear sigma model is given as 



H = j ^{^a(z) 2 + ^V0 Q (f) • V(p a (x) + \ [4>{x) 2 - v 2 ] 2 - h(j)o{x) \dx 



a=0 



(x) + V[<p(x)] \dx, (1) 



where 0(x) 2 = X)a=o rf'aix) 2 , and A, v, and h are constants. Note that ir a (x) is the conjugate field operator of 4> a (x) 
and should not be confused with the pion field operator. We determine the three model parameters, A, v, and h, so 
that they give the pion mass M T = 138 MeV, the sigma meson mass M a — 500 MeV, and the pion decay constant 
Jtc = 93 MeV in the one loop level in the "broken symmetry" vacuum state as in the previous study, i.e., A = 3.44, 
i; = 110 MeV, and /i = (103 MeV) 3 HHH. 

The squeezed states used as trial states have the following form: 



|$(<)> = expji^y {D a {x,t)4> a {g) -C a (x,i)Tr a (£))dx\ 

x A/-(;)expj^ J Mx)(\G { at\x^')~\Gab{S^'it)+i^ (2) 



.a, 6=0' 



Here |0) is the vacuum state, and G^aa(x,y) = {0\<p a (x)cj)a(y)\0) and G^{x,y) — for a ^ b. C a {x,t) and D a {x,t) are 
the c-number mean fields of (j) a (x) and n a (x), respectively. G a b(x, y, t) and n a fc(x, y, t) are the quantum correlation 
functions having isospin indices a and b (quantum fluctuation if both x — y and a — b hold), and the canonical 
conjugate variable for G a b(x, y, t), respectively. Thus, the trial state is specified by C a {x, t), D a (x, t), G a b(a?, y, i), and 
H a b(x,y,t). All of these variables are real. Af(t) is a normalization constant, whose explicit form is not needed in 
the calculation of the time evolution of observables. Its explicit expression is given in Ref. When there is no 

isospin-isospin correlation, the above squeezed state reduces to the direct product of the squeezed states with single 
isospin labels, a = ~ 3. The squeezed state include the coherent state as a special case. In addition, correlations 
can also be taken into account with the trial states. Thus, the trial states given by Eq. @ span a wide subspace in 
the physical Hilbert space and are expected to be able to describe a variety of quantum features of the system. 

The time evolution is determined by the time dependent variational principle: 



>Jm)\ift-Bm))dt = 0. (3) 
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The equations of motion in momentum space read 
C a (k,t) = -PG a (k,t)-MP(k,t), 



G ab (k,k',t) 
il ab (k,k',t) 
M<P$,t) 



= 2 



E (k\ [G aa ,(t)U a , b (t) + IW(*)Gv&(t)] \k'), 

a'=0 

3 3 

I E $\ G aUt)G-,\m') 2 E <*|na«'(*)n.'6W|A'> 

dldl 



a'=0 
,2 



a'=0 



i(27r) 3 fc 2 ,5 a ^ 3 (fc-fc') - -*'.*). 



m 2 C a (k,t) 



(2tt) 6 



4AC a (fc - f- ?, t) E C ^ t)C»(?. *) + 12AC„(fc - f- ?, ?, t) 



6=0 



-4AC a (£ -f-F,t) E t) + 8A E C b {k-l-?,t)G ia 0,t)) -hS aQ V 



6(#a) 



= ~™ 2 V+ I j^(l2X(C a (k-l,t)C a (l,t) + G aa {k-l,l,t)) 



= 8A 



4A E (C b (k-l,t)C b {l,t)+G bb (k-fXt))), 
dl 



(2tt) 3 



C a {k-l,t)G b {l,t) + G ab {k-l,l,t)), 

I — 



(4) 



where m 2 = 4Av 2 and V = Trl = J dx. C a (k,t) is 
the mean field for the 4> a field with momentum k, and 
G ab (k,k' , t) and n a ;,(fc, A;', £) are the correlation between 
modes with isospins a and 6, and momenta k and fc' 
(the quantum fluctuation for a — b and fc = k'), and 
the canonical conjugate variable for G ab (k,k' ,t), respec- 
tively. The one-point functions C a (k,t) and D a (k,t) 
and the two-point functions G ab (k, k' , t) and Il a ft(fc, k', t) 
are the momentum space representations of C a (x, t) , 
D a {x,t), G ab (x, x', i), and Ii ab {x, x', t), respectively. In 
Eq. Q), we have used the following notation, 



(k\H(t)I(t)\k') = 



1 



(2tt) 3 



H(k,k",t)I(k",k',t)dk". 



(5) 

Mf b \k,t) in Eq. @ originates from the nonlinear 
coupling term in the Hamiltonian and is given as follows: 

M^&ty^dx, 

3 



^(m) 



exp 



a'b'=0 



d 2 

* dz a dz b Z 



(6) 



2 a =C a (:r,*) 



where E7[z] represents the c- number potential that is ob- 
tained by replacing <ft a in the potential term of the model 
Hamiltonian V[<f>] by z a S R: 



U[z] = V[0] 



(a = 0-3). 



(7) 



The third equation in Eq. Q tells us that mode- 
mode correlations in momentum space arise through 

(2) ~* ~* — ♦ — * 

Ai al l(k — k ,t) with k k' even if there is initially no 
such correlation. Similarly, the isospin-isospin correla- 
tion arise due to the interaction term AIL (k — k , t) with 

In the numerical calculation, we assume that the sys- 
tem is in a cubic box with the volume V = L 3 and com- 
plies the periodic boundary condition, and that the lat- 
tice spacing in each direction is the same, d = -4, where 
N is the division number in each direction. Accordingly, 
the momentum fc, the square momentum k 2 , and the mo- 
mentum integral in Eq. Q are replaced as follows: 



k 
P 



2tt. 



dk 
(2^)3 



E4 . o [ITUi 
d^ Sm I AT 

i—x^ y,z 
, JV-1 JV-1 JV-1 



r 



E E E> 

n x — n.- n,— 



(8) 



with n 



B. Initial condition 

As the initial condition, we take the quench initial con- 
dition. In the quench scenario, the following assumption 
on the time evolution of the system is made: The mat- 
ter produced in relativistic heavy-ion collisions is initially 
thermalized and that chiral symmetry is restored. Then, 



4 



this matter goes through rapid cool down and the true 
vacuum spontaneously breaks chiral symmetry. In this 
process, however, the change of the effective potential 
is so rapid that the order parameter cannot follow the 
change and remains around the origin, where the mini- 
mum of the effective potential is in the chirally symmet- 
ric phase. P, 0,E3. Namely, in the quench scenario, 
the mean of the chiral fields as well as the correlation of 
the fields are the same as those before the sudden cool- 
ing (quenching), but only the sizes of the fluctuations of 
the mean fields and their conjugate fields are quenched. 
In order to realize the quench scenario in our numeri- 
cal calculation, we specify the initial condition for the 
mean fields and the quantum correlation and fluctuation 
as follows. 

For the mean fields, we have used the same prescription 
for the initial condition as in Refs. |ToL IT^| : at each grid 
site, the mean field variables for the chiral fields and their 
conjugate variables C a (x, 0) and D a (x, 0) are randomly 
distributed with the Gaussian weights with the following 
parameters, 



(C a (f,0)> 
(C a (x,0) 2 ) - (C a (x,0)f 
(D a (x,0)) 

(D a (x,0) 2 } - (D a (x,0)) 2 



0. 

0, 

d? ' 



(9) 



where <5 is the Gaussian width and D is the spatial di- 
mension. In relating the Gaussian widths of C a {x 1 0) and 
D a {x 1 0), we have taken advantage of the virial theorem. 
For detailed discussion on the use of virial theorem, see 
Ref. 0|- We use the Gaussian width S = 0.19v, which is 
the same value which has been taken in Ref. [10]. This 
value is so chosen that the fluctuation is small enough 
and the chiral symmetry is spontaneously broken. This 
is necessary to ensure that the quench scenario is simu- 
lated. This value of 5 is taken as a typical value and, of 
course, other choices are also possible. 

In order to fix the initial condition for the quantum 
correlation and fluctuation, we have assumed that each 
of the sigma and pion states in momentum space is in- 
dependently in a coherent state with a degenerate mass 
to - Accordingly, G ab (x, y, 0) and Tl ab (x, y, 0) are given 
as follows: 



G aa (x, y,0) 

Gab(^a) {X, y, 0) 

U ab (x,y,0) 



1 f A e^'C 5- ?) 



(2tt) 3 



dk. 



0, 
0, 



(10) 



with ujk = 



k 2 . We adopt m = 200 MeV. The 



two-point functions G a b(x, y, 0) and H a b(x, y, 0) are thus 
diagonal in both momentum and isospin spaces in the 
initial state. However, as previously explained, their 
off-diagonal components in both momentum and isospin 
spaces arise through the non-linear coupling term in the 



equation of motion Eq. Q , if the mean field is not trans- 
lationally invariant in coordinate space or symmetric in 
isospin space. We examine the effect of these off-diagonal 
components of the Green's function on the time evolution 
of the system in the next section. 



III. 



NUMERICAL RESULTS 



We evolve the mean fields and the two-point functions 
on a discrete lattice with the total length in the z direc- 
tion L — 64 fm and the lattice spacing d = 1.0 fm, which 
corresponds to the three dimensional momentum cutoff 
A = 1071 MeV. We assume translational invariance in 
the x and y directions, thus D = 1 in Eq. ©. We have 
confirmed that the total energy is conserved at least to 
an accuracy of less than 1 percent throughout the time 
evolution. 

To understand the role of the mode-mode and isospin- 
isospin correlations separately, we have performed the 
following four numerical simulations: 
Case i(a): The mode- mode and the isospin-isospin 
correlations are both included. 

Case i(b): The mode- mode correlation is ignored, but 

the isospin-isospin correlation is included. 

Case ii(a): The mode- mode correlation is included, but 

the isospin-isospin correlation is ignored. 

Case ii(b): The mode- mode and the isospin-isospin 

correlations are both ignored. 



A. Time evolution of the mean fields 

First we show the time evolution of the mean fields. 
Figure ^ illustrates the time evolution of the mean field 
of the third component of the pion field in the four cases. 
The value of the mean field is represented by different 
colors at each position and time. 

In case i(a) and case ii(a), where the mode-mode corre- 
lation is included, the formation of the domain structure 
of the pion field is clearly seen. Also it is observed that 
the size of the domains continues to grow beyond the 
time scale of the rolling down, ^a few fm, in these cases. 
On the contrary, the formation of only small domains is 
observed in case i(b) and case ii(b). In these cases, the 
mode-mode correlation is not included and short range 
fluctuation is dominant throughout the time evolution. 
No qualitative change in the behavior of the mean fields 
is observed after a few fm. These results imply that the 
mode-mode correlation is a key ingredient for the DCC 
domain formation. 

Furthermore, it is observed that the domain structure 
is larger in case i(a), where the mode-mode and isospin- 
isospin correlations are both included, compared with 
case ii(a), where only the mode-mode correlation is in- 
cluded. The appearance of the larger domains in case 
i(a) is due to the combined effect of mode-mode and 



5 




0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 

... , t [fm] 
n(a) 




0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 

t [fm] 



MeV] 




90.0 


1 


72.0 




54.0 




36.0 




18.0 




0.0 




-18.0 




-36.0 




-54.0 




-72.0 




-90.0 


MeV] 


-90.0 


172.0 




54.0 




36.0 




18.0 




0.0 




-18.0 




-36.0 




-54.0 


1-72.0 


1-90.0 




0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 
ii(b) ' [ fm ] 



■Si It 

0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 

t [fm] 




MeV] 


-90.0 


172.0 




54.0 




36.0 




18.0 




0.0 




-18.0 




-36.0 




-54.0 


1-72.0 


1-90.0 


MeV] 


-90.0 


|72.0 




54.0 




36.0 




18.0 




0.0 




-18.0 




-36.0 




-54.0 


1-72.0 


1-90.0 



FIG. 1: (Color online) Time evolution of the mean field of the third component of the pion field in cases i(a), i(b), ii(a), and 
ii(b). 
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FIG. 2: (Color online) Snap shots of the spatial correla- 
tions of the mean fields in case i(a) and case i(b). 



FIG. 3: (Color online) Snap shots of the spatial correla- 
tions of the mean fields in case ii(a) and case ii(b) . 



isospin-isospin correlations. Also it is observed that in 
case i(b), where the isospin-isospin correlation is included 
but the mode-mode correlation is not included, domains 
are slightly larger compared to those in case ii(b), where 
either correlation is not included. From these results, we 



find that the domain formation is most effective when 
both mode-mode and isospin-isospin correlations are in- 
cluded. 
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B. Spatial Correlation 

In this subsection, we show the spatial correlation func- 
tion C(r, t) defined by 



C(r,t) 



JC(x)-C(y)S(\x-y\-r)dxdy 
J\C(x)\\C(y)\S(\x-y\-r)dxdy 



(11) 



where C(x)-C(y) = £i=i Ci{x)d{y) and \C(x)\ = 



Tli=iCi{x) 2 ■ This is a measure for the closeness of 
the order parameters in isospin space at distance r in 
coordinate space. 

The snap shots of the spatial correlations at t = (the 
initial time), 5, 20, 40, and 60 fm are shown in Figs. |3 
and |3 for the four cases. The correlation functions are 
obtained by taking the ensemble average over 10 indepen- 
dent initial states. At the initial time t = 0, the spatial 
correlation length is almost equal to the lattice spacing 
in each case. This is so, because the initial mean fields 
are independently distributed at each grid site. However, 
difference in the correlation function appears as the time 
elapses. 

The correlation lengths in the upper figures are always 
longer than those in the lower figures in both Fig. and 
Fig. |31 This tells us the importance of the mode-mode 
correlation in the trial quantum states in the growth 
of the correlation length regardless of the presence of 
the isospin-isospin correlation. Moreover, the correla- 
tion length in case i(a) is longer than in case ii(a), while 
that in case i(b) is only a little longer than in case ii(b). 
The formation of larger DCC domains is induced by the 
combined effect of the mode-mode correlation and the 
isospin-isospin correlation. 

The duration of the correlation generation differs be- 
tween cases i(a) and i(b), and between cases ii(a) and 
ii(b). However, it is almost the same in cases i(a) and 
ii(a), and in cases i(b) and ii(b). In cases i(b) and 
ii(b), the formation of the correlation almost finishes by 
t = 5 fm, which is the typical time scale within which 
the rolling-down of the order parameter from the top of 
the effective potential is completed [l5j. We will explic- 
itly show this in the next subsection. The fact that the 
correlation increases beyond that time scale in cases i(a) 
and ii(a) strongly suggests that the mode- mode coupling 
is another driving force for the domain formation in addi- 
tion to the instability of the low momentum modes that 
exists during the rolling-down of the order parameter . 

As a possible mechanism of DCC formation, the para- 
metric resonance has been also proposed |l6j |. It is, how- 
ever, also understood within the framework of the mean 
field theory without correlations among modes; it is in 
essence a one-mode problem. Thus, the mode-mode cou- 
pling is yet another mechanism for DCC formation, and 
according to the result of our calculation, it has the most 
dominant effect on the generation of the field correlations. 
It is interesting to compare this result with the result of 
the classical simulation by Rajagopal and Wilczek p|. 
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FIG. 4: (Color online) Time evolution of the spatial average 
of the sigma field in cases i(a), i(b), ii(a), and ii(b). 



In their work, it was found that the amplification of low 
momentum modes lasts far beyond the typical time scale 
of the rolling-down of the order parameter. Their cal- 
culation is classical, but includes direct couplings among 
modes. This also exemplifies the importance of the direct 
mode-mode coupling in the domain formation. 



C. Time evolution of the sigma field 

In this subsection, we examine whether there is the 
collective oscillation of the sigma field, which is required 
for the DCC domain formation through the parametric 
resonance. For this purpose, we show the time evolution 
of the spatial average of the sigma field defined by, 



(Co(*)> = ~ 



Cq(x, t)dx. 



(12) 



Note that the zeroth component of the chiral field is the 
sigma field. 

Figure 0] shows the time evolution of the spatial aver- 
age of the sigma field in the four cases. This result was 
obtained by taking the ensemble average over 10 initial 
states. In each case the sigma field quickly rolls down 
from the top of the effective potential and approaches 
its ground state value within less than 5 fm. After that, 
the sigma field continues to oscillate as found in earlier 
studies |l|, l2| ■ The oscillation is damped more quickly in 
cases i(a) and ii(a), where the direct mode-mode coupling 
is taken into account. In case i(a), the oscillation almost 
disappears by t — 25 fm. On the other hand, according 
to Fig. [21 the correlation length of the pion field contin- 
ues to grow until about t = 40 fm. Thus, the collective 
oscillation of the sigma field is not necessarily needed for 
the DCC domain formation unlike in the parametric res- 
onance scenario. 
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FIG. 5: (Color online) Time evolution of the spatial average 
of the quantum fluctuation of the third component of the pion 
field, (G33(t)) in cases i(a), i(b), ii(a), and ii(b). 



FIG. 6: (Color online) Time evolution of the total particle 
numbers of the third component of the pion field N3 , its clas- 
sical part N§ 1 , and its quantum part JVJ in case ii(a). 



D. Quantum fluctuations 

In this subsection, we show the time evolution of the 
quantum fluctuation. 

Figure [5] depicts the time evolution of the spatial av- 
erage of the quantum fluctuation of the third component 
of the pion field, (G^t)), which is defined by, 

(G33W} = i J G 33 (x,x,t)dx. (13) 

In this figure, no event average is taken. However, since 
the spatial size of the system is large enough, the behavior 
of (Gn(i)) and {G22(t)) is expected to be similar to that 
of (G33(i)) due to the isospin symmetry. 

In the initial state, fluctuation is small. This is be- 
cause the initial state is chosen to be the direct product 
of the coherent states. The remarkable feature is that 
the quantum fluctuation is substantially amplified in the 
cases with the mode-mode correlation, i.e., cases i(a) and 
ii(a). On the other hand, such substantial amplification 
of the quantum fluctuation is absent in the cases without 
mode- mode correlation, i.e., cases i(b) and h(b). In ad- 
dition, it is noted that the duration of the amplification 
of the quantum fluctuation in cases i(a) and i(b) is of the 
same order as that of the correlation formation observed 
in Fig. [21 In the space of the trial states, the amplifi- 
cation of the quantum fluctuation implies squeezing of 
the states. Thus, the main part of the late time DCC 
domain formation occurs in concurrence with squeezing 
of the states. 

On the contrary, only little difference in the time evolu- 
tion of the quantum fluctuation is observed between cases 
i(a) and i(b), and also between cases ii(a) and ii(b) . This 
implies that the isospin-isospin correlation does not lead 
to further amplification of the quantum fluctuation. 

However, we have observed that there is a little differ- 
ence in the formation of correlations between cases i(a) 



and i(b), and also between cases ii(a) and ii(b), in Figs. 
12 and This suggests that there are at least two mecha- 
nisms that lead to the long range correlation of the pion 
fields: one that is accompanied by the squeezing of the 
pion fields and the other that is not. 



E. Pion particle numbers 

In this subsection, we calculate the particle numbers 
associated with the pion fields for cases ii(a) and ii(b). 

In order to analyze the particle numbers associated 
with the pion fields, we expand the field operator 4> a (x) 
and its conjugate operator Tr a {x) by the annihilation op- 
erator a a (k) and creation operator a^(fc) of the pion with 
the momentum k and isospin a, 



4>a(x) 



1 



(2tt) 3 

I 



I 



a a (ky k - s + ai(k)e- lk - s 
a a {ky % - 3 -a\{k)e- i% - d 



dk, 
dk, 



where tOk is the energy of the normal modes with the pion 



mass in the vacuum M„, u>k — y M% + k 2 . From these 
relations, we obtain 



a a (k) = 



e' ltcx dx, 
e itg dx. 



Then, the expectation value of the particle number op- 
erator of the pion with isospin a in the squeezed state is 
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FIG. 7: (Color online) Time evolution of the total particle 
numbers of the third component of pion field N3, its classical 
part N§ 1 , and its quantum part N% in case ii(b). 



given as follows: 



n a (k,t) 
nl(k,t) 



mt)\ai(k)a a (k)\m) 
n c a l (k,t)+nl(k,t), 

^C* a (k, t)C a (k, t) + ^j- k D*(k, t)D a (k, t) 

l -(C;(k,t)D a {k,t)-C a (k,t)Dl(k,t)), 
~~^~G aa (k, k, t) 

^-(iG- Q 1 (fc,fc,i) + 4(fc|n aa G aa n aa |fc)) 



+ i(k\(G aa U aa -U aa G aa )\k), 



(14) 



where we have divided n a (k,t) into the classical part 
n^ l (k,t) and the quantum part n^(k,t); n^ l (k,t) in- 
cludes the mean field and its time derivative, C a (k, t) and 
D a (k,t), and n*(fc, t) includes the quantities associated 
with the quantum fluctuation. On the right hand side 
of Eq. I|14|) . n o sum over a is taken. In particular, no 
off-diagonal quantity with regard to the isospin indices 
appear on the right hand side of the last equation of ifh4"|) 
because no isospin-isospin correlation exists in either case 
ii(a) or case ii(b). In addition, we define the total parti- 
cle numbers as the sum of the particle numbers over all 
modes, 

N a(t) = 



k 

K(t) = 5>«(M). 



(15) 




We have calculated these particle numbers by taking 
the ensemble average over 10 different initial field config- 



FIG. 8: (Color online) Snap shots of the two-point function 
in momentum space, IG33 (k, k' , in case i(a). Note that the 
scale of the color legend at the initial time t — is different 
from the others. 



urations. Figs. El and show the time evolution of Na(t), 
N^ l (t), and iVg(i) in cases ii(a) and ii(b), respectively. 

We observe that the total particle number Ns(t) 
rapidly increases within the typical time scale of the 
rolling-down of the order parameters (a few fm) in both 
cases. At first, the classical part N^ l (t) is dominant in 
both cases. In case ii(a), where the mode- mode correla- 
tion is included, the classical part Nf(t) begins to de- 
crease and the quantum part A|(i) begins to increase 
after a few fm, while the sum remains almost unchanged. 
The typical time scale of the increase of the quantum 
part is again close to that of the formation of the long 
range correlation. This clear separation of the stages in 
the time evolution of N§ l (t) and N!f(t) also shows that 
the main mechanism for the domain formation is not the 
rolling-down of the sigma field and the associated insta- 
bility of pion fields with low momenta. On the other 
hand, in case h(b), where less correlation results, the 
classical part of the total particle number N^ 1 (t) is dom- 
inant throughout the time evolution. 



F. Off-diagonal components of the two-point 
functions in momentum space and isospin space 

In this subsection, we show the time evolution of the 
off-diagonal components of the two-point functions in 
momentum space and also in isospin space, which corre- 
spond to the mode-mode correlation and isospin-isospin 
correlation, respectively. 

Figure El shows the absolute value of the two-point 
function within the third component of the pion field 
in momentum space, \Gas(k, k', t)\, in case i(a). In Fig. 
El no event average is taken. As explained in Section 
ITT1 in general cases where there is no translational in- 
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FIG. 9: (Color online) Snap shots of the two-point function 
in isospin space, (|G a &(i)|), in case i(a). Note that the scale 
of the color legend at the initial time t — is different from 
the others. 

variance, the off-diagonal components of the two-point 
functions develop and grow as the system evolves in time 
even if they are initially absent. According to Fig- El the 
off-diagonal components indeed appear as time elapses. 
Although the diagonal components are always dominant, 
the off-diagonal components are not negligible, in partic- 
ular, at low momenta. This existence of the off-diagonal 
components is crucial for the DCC formation and the 
evolution of the quantum fluctuation, as we have already 
demonstrated. 

FigureElrepresents the time evolution of the two-point 
functions in isospin space in case i(a). In this figure, the 
spatial average of their absolute value, 

(\G a b(t)\) = i J \G ab (x,x,t)\dx, (16) 

is shown. No event average is taken in Fig. either. 
At the initial time, the two-point function is diagonal in 
isospin space as in momentum space. The off-diagonal 
components of the two-point function emerge and de- 
velop as time elapses. The absolute value of the off- 
diagonal components evolves up to about one fourth of 
that of the diagonal components by t ~ 25 fm, and then 
their relative strength starts to decrease. 



G. Dissipation of the sigma field 

Finally, we explore whether and how the dissipation of 
the sigma field is included in the approximation. For this 
purpose, we chose the following initial condition. The 
mean field of the sigma field Co was randomly distributed 
according to the same Gaussian form as Eq. © except 
that the first equation was replaced by (Co) = 60 MeV. 
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FIG. 10: (Color online) Time evolution of the spatial average 
of the sigma field. 

The pion fields and their conjugate fields were set to zero 
at each lattice point. As for the quantum fluctuation 
and correlation, the same initial condition as in Section 
III Bl i.e., those of independent coherent states with the 
degenerate mass mo = 200 MeV, were used. Ensemble 
average was taken over 10 initial configurations as before. 
In the calculation, both mode-mode and isospin-isospin 
correlations are taken into account, i.e., case i(a). 

Figure ^] represents the time evolution of the spatial 
average of the mean field part of the sigma field. The 
result tells us that the oscillation of the sigma field is 
damped quickly and it approaches a constant. Since ini- 
tially the pion fields and their conjugate variables are 
all zero, the mean field part of the pion fields remains 
zero throughout the time evolution. On the contrary, 
the quantum fluctuation of the pion fields is not zero 
initially and, moreover, it is amplified substantially as 
shown for the third component in Fig. 1111 while that 
of the sigma field changes only little. This amplifica- 
tion does not exist in the translationally invariant case. 
Physically, the damping of the oscillation of the sigma 
field and the growth of the quantum fluctuation of the 
pion fields correspond to the decays of the sigma into two 
pions. As observed from Figs. 1101 and ITT1 the growth of 
the fluctuation of the pion fields stops approximately at 
the same time as the oscillation of the sigma fields ends. 



IV. CONCLUSION 

In this paper, we have studied the dynamics of chi- 
ral phase transition in spatially inhomogeneous systems 
with mode-mode correlations in the framework of TDVA 
with squeezed states for the 1+1 dimensional geometry. 
Quantum effects and spatial inhomogeneity are both in- 
corporated. 

In the case with the mode-mode correlation, the DCC 
domain formation continues beyond the time scale of the 
rolling-down of the order parameter. The quantum flue- 
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correlation on the domain formation of DCC for the 
first time. We found that the isospin-isospin correlation 
makes the DCC domain structure larger, but its effect is 
less conspicuous than that of the mode-mode correlation. 
The isospin-isospin correlation has almost negligible ef- 
fects on the time evolution of the quantum fluctuation. 
Thus, the isospin-isospin correlation does not lead to fur- 
ther squeezing of the states. 

We plan to include realistic geometry and dynamics 
such as expansion, and carry out calculations in 2+1 or 
3+1 dimension in future work. They are indispensable 
for the quantum mechanical understanding of the DCC 
formation in ultra- relativistic heavy ion collisions. 



FIG. 11: (Color online) Time evolution of the quantum fluc- 
tuation of the sigma field (00) and the third component of the 
pion field (33). 

tuation is also amplified with approximately the same 
time scale. This implies that squeezing of the states 
takes place. On the other hand, no significant DCC do- 
main formation or squeezing of the states was observed in 
the case without the mode-mode correlation. Thus, it is 
concluded that the mode-mode correlation plays the key 
role in the formation of DCC domain. This mode-mode 
correlation has not been taken into account in preceding 
quantum calculations for DCC formation. In addition, it 
was pointed out that this mechanism is different from the 
previously proposed parametric resonance amplification. 
We have also examined the effect of the isospin-isospin 
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